a density peak are changed when small scale substructure induced by dynamical friction 
are taken into account. The peak mass obtained is compared to the result of Peacock 
& Heavens (1990) and to the peak mass when dynamical friction is absent to show how 
these effects conspire to reduce the mass accreted by the peak. 
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1. Introduction 

The origin and evolution of large scale structure is today the outstanding problem 
in Cosmology. In this context, a remarkable role has been played by the standard Cold 
Dark Matter (CDM) model. 

The main features of the model are a scale invariant spectrum of density perturbations, 
P(k), of the type Harryson-Zeldovich, growing under gravitational instability and the 
supposition that the bulk of the material in the Universe should be non baryonic but un- 
der the form of Weakly Interacting Massive Particles (WIMPS) (Efstathiou 1990, Kolb 
& Turner 1990). The only addition to the model, required by observation, is the question 
of biasing, i.e. the fact that the distribution of group and clusters of galaxies is different 

J- — from the distribution of cosmic material on large scale, the first being more clustered 
than the second. The physical origin of such a biasing is not yet totally clear even though 
several mechanisms have been proposed (Rees 1985, Dekel & Rees 1987, Colafrancesco 

^ et al.l995). 

*^ In biased galaxy formation theory, it is assumed that structures of size Rf form around 
I the local maxima of the density field, smoothed on the filtering scale Rf. The linear 
^ density perturbations evolve towards the non linear regime because of gravitational in- 
stability and collapse when the average density is, (5 ~ 1. 

A peak of the primordial density fiuctuation field has a probability to form a protostruc- 
ture proportional to its central height: 

i (1) 

where (5(r) is the overdensity in the radius r and ao(Rf) the mass variance windowed by 
a filter of scale Rj. It is assumed that only peaks with > i/^ = = ^^^^^ -) form 

structures. 

The assignement of a mass, Mpk, to a peak of the filtered density perturbation field, is a 



Heavens 1990; Lucchin & Matarrese 1988; Ryden 1988; Colafrancesco et al. 1989) have 
given different solutions. 

The first and simplest is to assign to the peak the mass obtained from the volume of the 
filter function, Rf: 

Mpk = (27r)3/Vci?/ (2) 

(Efstathiou & Rees 1988) where pc is the critical density. This means that the mass 
contained within a given volume is produced by the filtering process and we suppose 
that an object of lenght Rf collapses when the mean overdensity in Rf is, as previously 
written, (5^ ~ 1. 

If otherwise we use a spherical filter, the mass is given by: 

Mpk = y^/Pc (3) 

An alternative to those given is to estimate Mpk by modelling a peak as a triaxial ellipsoid 
(Peacock & Heavens 1990; Colafrancesco et al. 1989) obtaining a distribution of mass for 
a given i>: 
Ait 

Mpk = -YPcaia2a3 (4) 

where ai, 02, 03 are the semiaxes of the ellipsoid. In terms of v, i>c, the parameter of 
ellipticity e, the parameter of oblateness p, the local curvature parameter of the peak and 



the mass M^, = ^R^Pc, where R^, cx Rf we have: 



Mpk{:^, V,, M; e,p, x) = (1 - fpu -fl^—fl^ { [(1 + pf - <de^] (1 - 2p)x^}-'''' (5) 

Vc 1 

(Colafrancesco et al. 1989), where /p^ is a term simulating the mass loss by the peak due 
to tidal interactions with neighbours, 7 is a parameter that can be expressed in terms of 
the spectral moments (see Bardeen et al 1986): 



,2 



k^dk 



P(k,t)k^^ (6) 



27r2 

Integrating Eq. (5) over e and p we have: 

Mpk=il-fpk^] pij,iyj)M, 0<7i.<1.5 (7) 

where p(j,ji>) is a function given by Colafrancesco et al.(1989) (Eq. 6). 
Peacock & Heavens (1990) give another expression for the peak mass: 

23/2(47r)pei?g ^ns (fi\ 

= ^3 + (M)3/2 0.5 < 7 < 0.8 (8) 



cess of the density field and on the filtering scale, which in turn depends on the physical 
process which set the threshold i>c. The exact analytic form of the filtering function is 
not known and it is often taken to be a Gaussian. At the same time there is not a precise 
relation between the filtering scale, Rf, and the mass of a class of objects (as we have seen 
in the Eq. (2), Eq. (3) and following). A more complex method, trying to overcome these 
problems, is that of Appel & Jones (1990). It is based on the observation that, during 
the process of filtering of the density field when Rf is increased, some regions lying over 
the selected threshold grow in size while others decrease. The ones that decrease in size 
will vanish for a certain value of Rf. The point of vanishing is the position of the local 
maximum and Rf is its radius. In this way the mass of a peak is chosen by the same 
procedure of filtering. 

Finally, there is another procedure to define a peak mass and it is connected to the def- 
inition of a binding radius for the structure (Hoffmann & Shaham 1985). It is directly 
connected to the process of collapse of a density perturbation. 

As shown by Antonuccio-Delogu & Colafrancesco (1994) the evolution of structures like 
Clusters of galaxies is infiuenced by the dynamical friction induced by the substructure 
present in them in accordance with CDM model. Consequently the binding radius and 
the mass of a density peak must be infiuenced by dynamical friction. These observations 
lead us to recalculate the peak mass taking into account this dissipative effect. 
In this paper we use the definition of mass peak given by Hoffmann & Shaham (1985) 
to calculate the peak mass and we show how it must be changed due to the presence 
of substructure in protostructure, inducing dynamical friction. Finally, we compare this 
result to that given by Peacock & Heavens (1990). 

2. Mass peak from binding radius 

There are at least two ways to assign a radius to a peak of density: a statistical one, 
(Ryden 1988), and a dynamical one, (Hoffmann & Shaham 1985). 

In the first criterion the binding radius of the region, rj, is given by the solution of the 
equation: 

< 6{r) >=< (6{r)- < 6{r) > f >^l'' (9) 

At radius r << rj the motion of particles is predominantly toward the peak while when 
r >> rj particles are not bound to the peak. 



shell collapses in a time, Tco, smaller than the age of the universe to'. 



Tco{r) < to 



(10) 



This last criterion, differently from the previous one, contains a prescription particularly 
connected with the physics of the collapse process of a shell. It supposes that only the 
regions that have collapse time less than a Hubble time form structure. The time of 
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Fig. 1. Variation of the binding radius rt with j/. The solid line is the binding radius in absence 
of dynamical friction, while the dashed line is the same as in presence of dynamical friction. The 
filtering radius used is , i?/ = 1 Mpc 



collapse, Tco(r), can be obtained using the equation: 

TeoW=^(l±I^ (11) 

Gunn k Gott (1972). 

To obtain the average density b we used the equation given by Bardeen et al. (1986): 



8{r) = A' 



vi{r) 



.^(0)1/2 ^^(0)1/2(1 _ ^2) _ 



(12) 



^(r) is the correlation function of two points, 7 and i?* two constants obtainable from the 
spectrum (see Bardeen et al. 1986) and finally 0(ji>, 7) is a function given in the quoted 
paper (Eq. 6.14). The average density, 6, inside the radius r in a spherical perturbation 
is given by: 



We calculated the time of collapse, Tco(r), using Eq. (11) with the density profile given 
in Eq. (12). We repeated the calculation of Tco(r) for 1.5 < < 3 and we applied the 
condition given in Eq. 10 to the curves Tco(r) previously obtained. The result is the plot 
in Fig. 1 for the binding radius rj versus i>. Using the plot ri, — i> and the equation: 



finally we obtained the peak mass. As expected, it is an increasing function of i>. 
3. Mass of a peak in presence of dynamical friction 

The calculation of the peak mass, given in the Sect. 2, does not take into account the 
substructure present into a collapsed perturbation. In presence of it the collapse of a 
structure is modified (Antonuccio-Delogu & Colafrancesco 1994) and, as a consequence, 
the mass accreted by the structure is modified too. 
The time of collapse in presence of dynamical friction is given by: 



where Aq is a function of the number of peaks and of its average mass (Eq. 23 in 
Antonuccio-Delogu & Colafrancesco 1994), while c(6) is a function of the overdensity 
6 (Eq. 27 in the quoted paper). 

We calculated the collapse time Tc(r) for a given value of using Eq. (15) and the density 
profile Eq. (12). We repeated the calculation procedure made in the previous section to 
determine the binding radius, rj for 1.5 < < 3. Using the plot ri,-i> and the Eq. (14) 
we finally obtained the peak mass, Mpk in presence of dynamical friction. In Fig. 2 we 
compared this result with the peak mass in absence of dynamical friction and with the 
equation of Peacock & Heavens (1990). 

The peak mass in absence of dynamical friction is comparable to that given by Peacock 
& Heavens (1990) while, as expected, dynamical friction reduces the peak mass with 
respect to systems in which it is absent. The definition of the peak mass is infiuenced by 
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Fig. 2. Variation of tlie pealcs mass Mph witli v. Tlie solid line representes the curve derived 
by Peacok & Heavens given in text; the dashed line rappresentes Mph in absence of dynamical 
friction and the dotted line rappresentes Mph in presence of dynamical friction 

the presence of substructure and then from dynamical friction. The total mass density 
associated with peaks is consequently changed with also the construction of the mass 
function following the Press & Schechter (1974) prescription. 

This analysis shows that dynamical friction reduces the mass bound by a peak with re- 
spect to that they should acquire if this dissipative effect was not present. 
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